Concept of local polaritons and optical properties of mixed polar crystals* 



Lev I. Deych^ Alexey Yamilov,-'- and Alexander A. Lisyansky^ 
^Department of Physics, Seton Hall University, 400 South Orange Ave, 
South Orange, NJ 07079 

. . . ^Department of Physics, Queens College of CUNY, Flushing, NY 11367 

O I (February 1, 2008) 

, The concept of local polaritons is used to describe optical properties of mixed crystals in the 

' frequency region of their restrahlen band. It is shown that this concept allows for a physically 

^ ' transparent explanation of the presence of weak features in the spectra of so called one-mode crystals, 

^ , and for one-two mode behavior. The previous models were able to explain these features only with 

the use of many fitting parameters. We show that under certain conditions new impurity- induced 
polariton modes may arise within the restrahlen of the host crystals, and study their dispersion laws 
and density of states. Particularly, we find that the group velocity of these excitations is proportional 
to the concentration of the impurities and can be thousands of times smaller then the speed of light 
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I. INTRODUCTION 



The optical properties of mixed polar crystals have been attracting a great deal of interest since the 1950's. Main 
\^ i efforts have been directed to experimental and theoretical studies of the concentration dependencies of fundamental 
transverse (TO) and longitudinal (LO) phonon modes, and to the fine structure of spectra in the frequency region 
Q . between them {Restrahlen band). Reviews of earlier experimental and theoretical works in this area can be found 
Q in Refs. 0J|. In spite of the disordered nature of mixed crystals, it is usually possible to observe both TOpaad LO 
modes of pure crystals at both ends of the concentration range, as well as features associated with disorder .ETeI This 
' is usually done in reflectance or transmittance experiments by identifying the maxima of Im(e) and — /m(l/e) with 
^ TO and LO modes respectively, where e is the dielectric function.cl 

^vq ' The term mixed crystals is usually applied to materials in which the concentration of each component is large, so 
I that no component can be considered as a system of independent impurities. At smaller concentrations the dynamics 
• of impure systems is described in terms of defectj-^tates, which are either localized or quasilocalized depending on 
\l ' whether they fall into forbidden or allowed bandaJ B of the pure system. The transition from low to high concentration 
behavior occurs when the average distance between defects becomes of the same order as the localization length of 
the individual states. In the case of the local phonon states, the localization length is of the order of a few lattice 
constants, and the respective transition concentration is of the order of ~ 10^^ cm^^. Most of the studies of mixed 
■ crystals were focused upon combinations of alkali halides, II/IV and III/V group polar binary crystals ABi-pCp, 
] where p is concentration. The B atoms are substituted with the atoms C from the same column of tiieLPeriodic Table. 
I , It was found that there were several patterns that TO and LO modes follow in mixed crystals.ETLrcl'Ej In the type I 
' mixed crystals, also called one-mode crystals, the frequencies of the TO and LO modes evolve smoothly and almost 
^ , linearly with the composition parameter p from their values in AB to those in AC. In some one-mode systems (e.g., 
O ■ Ba/SrF2, Ca/ SrFj^3}3) an additional pair of weak modes inside the absorption band is observed. 

Type II or two-mode mixed crystals exhibit qualitatively different behavior. Two distinct restrahlen bands are 
usually observed in the crystals of this category. For arbitrary concentration, the width of each of these bands, 
^Ilo ~ '^TOi is approximately proportional to the composition parameter p for the dopant related band, and to 1 — p 
for that of the original material. Two-mode behavior is characteristic for the crystals formed by elements from III 
and V groups of the Periodic Table (e.g. GaAs/P and the like). l— . 

AnotheMnode beitavipr could not be fit experimentally either to one- or to two-mode patterns (e.g. Ga/ InP)^, 
PbSe/Te^ K/i^WPta etc.). This group of mixed crystals was called Type I-II or one-two-mode systems. In these 
systems, one mode behaves as in Type I, while the other behaves as in Type II, which is why Type I-II is also called 
mixed mode. 
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To explain different types of behavior in mixed polar crystals, a number of theoretical models have been proposed. 
In the-JEr, experiments, modes with wave numbers g are excited. The Random Element Isodispacement model 
(REI)ll3H tajkes advantage of this fact, assuming that the sublattices vibrate in phase with q = 0. With some 
modificationsO after including the local field, and assuming a dependence of the force constants upon the composition, 
the modified REI model (MREI) was successfully used to fit two-mode crystals but encountered some problems 
in the one-mode materials. The fine structure of some mixed crystals was reproduced along the lines of REI by 
means of a considerable increase in the number of fitting parameters (cluster modelEa). A different approach made 
use of the Green's function formalism in order to treat vibrations in the mixed crystals; averaging over disorder was 
performed either in a sinmJe virtual crystal approximation (VGA) or with the use of the much more elaborate coherent 
potential approximatiorO (GPA). Interaction with the electromagnetic field in these approaches was included after 
the averaging procedure was performed. The mean field approximation treats the crystal as if it has atoms of the 
same sort with some average properties. It reproduces well the one-mode behavior, but could not possibly explain 
the two-mode systems. GPA combined with an electrostatic treatment of the electric field was used to fit a broader 
range of experimental data (see, for instance, Ref. |20| ). 

An important theoretical problem many researchers have focused upon is to find a simple criterion for different 
mode behavior, based upon dynamical properties of the crystals. LucovskyO originally considered the significance of 
whethecjthe original absorption bands overlap or not, which turned out to be a rather rough criterion. Harada and 
Naritall3 proposed a criterion based on relations between MREI constants (see also review papers in Refs. It 
turned out that one-mode and two-mode crystals can be described using quite simple models with only a few fitting 
parameter, such as VGA and MREI, respectively. Lucovsky et al. in Ref. ^ made an attempt to establish a connection 
between the behavior of a mixed crystal and the existence of the localized modes at small values of concentration 
parameters p ov 1 — p. It was suggested that the two-mode behavior is associated with the existence of local impurity 
states at both ends of the concentration range, while the one-mode regime occurs when such local states do not arise. 
These ideas, however, could not be applied to the mixed-mode behavior since the local phonon stateSpiarise either 
in the gap between acoustic and optic branches or above the optic branch of the phonon spectrum.ETQ Therefore, 
the mixed-mode crystals required much more elaborate models, where agreement with experimental data can only 
be achieved by increasing the number of fitting parameters. The same was true for weak features of spectra of the 
one-mode crystals. 

One of the goals of the present paper is to put forward a simple physical picture explaining the mixed-mode behavior 
based upon the concept of local polaritons and an impurity-induced polariton band. The local polariton is a state that 
arises due to an impurity with frequencies inside the restrahlen band of a crystal. These states have an electromagnetic 
component which is localized around the defect, because it cannot propagate in the region with a negative dielectric 
parameter. A local polariton state within the restrahlen region was first observed experimentally in Ref. 21 as a 



weak feature in reflectance and Raman spectra of GaAs doped with Sn, Te or S. These states were found to be 
associated with the local change in susceptibility due to localization of an electron around the dopant. An interaction 
between the localized electron and LO vibrations of the crystal gives rise to local LO phonons, therefore these local 
states have a complicated structure that involves interactions between electrons, phonons and-electromagnetic field. 
Theoretically these states, however, were only considered in the electrostatic approximation,Eil and therefore these 
early observations did not lead to the concept of local polaritons, which requires that the retardation be taken into 
account. 

A new wave of interest in the optical properties of impure crystals in the restrahlen region occured recently and 
was due to a general interest in systems with depleted, or altered electropaagnetic density of states (DOS). The 
primary examples of such systems are photonic crystals,Ej and microcavitieso The restrahlen region in polar crystals 



was considered from this new perspective independently in Refs. 24 -2^, where the concept of local polaritons was 
introduced. The local state considered in Ref. ^ arises due to an impurity atom with optically active electronic 
transitions, which interacts with host atoms through the electromagnetic field only. With the transition frequency 
inside the restrahlen band, this atom fornis, a local atom-radiation bound state. To some extent these states are 
similar to the ones observed by Dean et al.^ though an interaction with phonons was left out in Ref. |2j. A different 
type of local polariton was considered in Refs. where it was shown that regular isotopic impurities can give rise 

to local states with frequencies inside the restrahlen region. The interaction with the retarded electromagnetic field 
is responsible for the localized electromagnetic component of the states. The detailed analysis of the local polariton 
mode in a three dimensional sodium-chloride-like structure (BGG) crystal in the case of diagonal and off-diagonal 
disorder was carried out in Ref. ^ It was proven that, because of retardation effects, a local polariton mode splits 
off the bottom of the TO-LO gap for an arbitrarily small strength of the defects even in three dimensional systems. 

When considering local polaritons, the authors of Refs. p5|-p7| assumed that the restrahlen is a spectral gap not 
only for electromagnetic excitations but for phonons as well. This is indeed the case for some polar crystals. It occurs 
more frequently, however, that the restrahlen region is devoid of transverse optic phonons, but is still filled with LO 
modes. In this case the interaction between local polaritons and LO phonons would make the former quasi-stationary. 
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Whether the local polaritons survive this interaction depends upon their life-time, and we shall address this question 
in the present paper. At the same time it is useful to note that the states observed in Ref. |2l| reside in the frequency 
region where the density of LO modes is especially large. The fact that they remained observable allows for optimism 
that the local polariton states of Refs. ^-p7| could also survive provided that DOS of LO is not too large. If this is 
the case, then local polaritons due to isotopic impurities can be invoked to explain the mixed- mode behavior of mixed 
crystals and weak features in the one-mode systems without having to introduce tens of fitting parameters. 

These ideas can, for example, be used to discuss the one-two-mode mixed crystal Gao.7oIno,3oP. Pure GaP has 
a complete polariton gap in all directions, and experiments of Ref. ^ clearly demonstrated the polariton band 
associated with In atoms inside iiie restrahl of GaP even at room temperature. Polariton branches inside restrahlen 
were also observed in pure CuCLEj This material demonstrates a very peculiar behavin&Jsecause Cu atoms can occupy 
several non-equivalent positions, thereby producing defects even in a pure m£itp:ial.Q'c3 These defects were shown to 
be responsible for a new TO mode with a frequency in the main restrahlen!u^ These experiments were originally 
described using a phenomenological two-oscillator model, but they perfectly fit to the idea of local polariton sates. 

Local polaritons, however, can take us further than an explanation of the old reflectivity experiments. When 
retardation is taken into account, the local polaritons give rise to a new transmitting channel for electromagnetic 
excitations, via impurity-induced polariton band. Experimentally, such a band can be observed in transmission 
steady-wave and time-resolved experiments. The majority of the old experiments mentioned above dealt with reflection 
spectra. Dips in the reflection coefRcient in those experiments were associated with impurity-induced absorption, while 
they could have actually occured because of enhanced transmission via the impurity-induced polariton— band. The 
effect of the enhanced transmission inside the restrahlen region of a polar crystal was observed in CmCLEJ 

We showed analytically that lof^ polaritons give rise to a strong resonance enhancement of the transmittance at 
the frequency of the local mode,y'E3 using a one-dimensional chain of dipoles. In the system without dissipation, the 
maximum value of the transmission coefficient reaches unity, although peak width depends exponentially upon the size 
of the system. Numerical simulations for a finite conceiitration of impurities showed that a localized mode develops 
into a conduction band inside the original polariton gap.EHl The growth of the local state into a conduction band in the 
system of optically active defects was also demonstrated in Ref. An analytical treatment of the impurity-induced 
polariton band in the dispersionless one dimensional dipole chain, with_cpncentration of defects varying in the range 
< p < 1 was proposed in Ref. |3j. Using the microcanonical methodcS we were able to calculate DOS due to the 
defect subsystem and the extinction (localization) length. We found that the impurity-induced band has a number of 
interesting properties. The group velocity of electromagnetic excitations, propagating in this band, for instance, was 
found to be proportional to the concentration of the impurities, and can be significantly smaller than the speed of light 
in vacuum. The results of analytical calculations were in agreement with numerical simulation in the presence of the 
spatial dispersion. The knowledge obtained from analytical dependencies of DOS in the one-dimensional model gives 
us insight into a treatment of three dimensional mixed polar crystals. The spatial size of local polariton states may 
be as large as the EM wavelength, meaning that even at a very low impurity concentration, ^ 10^^ cto~^, the local 
polaritons significantly overlap. This fact allows us to develop a continuous approximation for calculating properties 
of the polariton-induced band. Applying this approximation to the one-dimensional case, we can use our previous 
one-dimensional results to test the suggested approach. 

The structure of our paper is as follows. In Section II, the continuous approximation for the mixed crystals, based 
upon the concept of local polaritons, is introduced. In Section III, we estimate the life-time of the polariton local 
states due to interaction with LO modes. We show that under realistic assumptions about the density of states of 
LO modes, this life-time is compatible with that caused by inelastic relaxation, and may be long enough for local 
polaritons to survive. Section IV deals with the properties of the impurity-induced polariton band (dispersion laws 
and density of states of the respective excitations are studied). We also consider in this section the concentration 
dependence of poles and zeroes of the dielectric function, and a simple explanation of the mixed-mode and one-mode 
(with weak structure) behaviors is proposed. Section V is devoted to investigation of the reflection and transmission 
spectra of the semi-infinite mixed crystals and a slab of finite width. In Section VI we calculate scattering length 
of the impurity induced polaritons due to the fluctuations of concentration, and demonstrate why this scattering is 
significantly less important for impurity-induced polaritons then for regular polariton branches. This section provides 
an additional justification for the main results of the paper. 

II. MODEL 



In our recent work, Ref. 34, we considered development-of an impurity-induced polariton band in a one-dimensional 
model, and used the so-called microcanonical methodtJ in order to analytically calculate the complex Lyapunov 
exponent of the system. The approximation used turned out to be in a very good agreement with numerical results, 
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and was shown to be equivalent to a continuous medium approximation. The latter can easily be generalized for 3-d 
systems, and below we develop an approach to the impurity-induced polariton band in 3D mixed polar crystal using 
the following fundamental assumptions. 

In order to model optic vibrations of the system under study, we introduce two subsystems of oscillators with 
different frequencies. One subsystem represents an optic mode of the host atoms. The interaction between this mode 
and light leads to the polariton gap between TO and LO frequencies of the pure crystal. The second subsystem 
introduces vibrations of impurities, and it is assumed that its frequency belongs to the restrahlen of the host crystal. 
This is a crucial assumption of the model, since interaction between light and this mode brings about local polariton 
modes. The significance of this assumption rests upon two ideas. The first one goes back to Lucovsky,tl who was the 
first to connect modes of mixed crystals with local impurity modes. The second one is the concept of local polaritons, 
which are local states with fpjouencies in the restrahlen of the host. In the one-dimensional case we have confirmed 
analytically and numericallytHIEil that the presence of the local polaritons gives rise to an impurity polariton branch, 
when the concentration of impurities grows. 

The frequencies of both oscillators are, in general, complex- valued. The imaginary part of host vibrations is due 
to anharmonicity, while the impurity mode has two sources of decay. Besides usual relaxation, local polaritons can 
acquire an additional imaginary part due to interaction with LO phonons which in some cases can fill the restrahlen. 
We will consider the contribution of this interaction into the life-time of local polaritons in the next section, and the 
general effects of relaxatio n upon reflection spectra of mixed crystals will be discussed in Section V. 

It was found in Refs. p5|- ^7[|34| that for all frequencies outside of the immediate vicinity of the TO boundary of the 
restrahlen, the localization length Zq of local polaritons is of the order of magnitude of the wavelength of the incident 
light ~ 10~^ cm. Even for residual concentration n ~ 10^^ cm~^ of the impurities, nl^ is large. The individual 
states significantly overlap and a macroscopic volume containing many impurities can still be much smaller than the 
localization length. This fact allows us to develop a continuous medium approximation similar to the usual one used 
to treat long-wave excitations£3 but this time we use it to treat the sub-system of impurities. We would like to 
emphasize again that this consideration is reasonable for the polariton band only because of the long-range nature of 
the localized electromagnetic component of the local polaritons. 

The microscopical Hamiltonian describing optic modes of our system is: 

r ^ ^ r' r,r' 



where /^o, Mi ^-nd Ur, Vr' are the reduced masses and relative displacements of the ions in pure AB and AC crystals, 
respectively. We will assume that the presence of atoms C does not affect the order in the A sublattice. The summation 
runs over spatial index r for host atoms and over r' for the impurities. The third term in the above expression is the 
potential energy. In the harmonic approximation it can be written as 

$ (ur, VrO = <^ (r - n) u>i:^ + <^^^} (r' - r'l) 

+ $(oi)(r-r')<f^,. (2) 

In order to avoid unnecessary complications, we assume no spatial dispersion in our model Hamiltonian and 
isotropy of force constants 

$(Ur,VrO =*OU?+$lV?,, (3) 

where $0 and 'I'l are force constants describing the interaction of the ions in the AB and AC lattices, respectively. 
As follows from the form of expression the contributions due to host and impurity ions in the Hamiltonian (|l|) 
can be separated 



r 



1 / dUr \ ^ 1 ^ 2 

1 f dw,,V 1 



(4) 



We have several parameters of the dimension of length in our system: a lattice constant, a; an average distance 
between defects, I ^ n~^/^, which depends on concentration; and the localization length and the wavelength of the 
incident light, which are of the same order, Iq ^ \ ^ 10^^ m. We shall assume that the concentration is such that 
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n ^ Iq^. When this condition is satisfied, the individual defect states overlap, and the dynamical properties, apart 
from those at the fluctuation band-edge, do not depend significantly upon the particular arrangement of the defects 
in the crystal. Hence, one can introduce a smoothing parameter Isv such that I <^ Isv ^ ^o- The macroscopic 
volume associated with this parameter is such that it contains macroscopically many impurities, but there are still 
no significant changes in spatial distribution of vibrations over this length: Ur+R, Vr+R ~ Ur, Vr for R ^ Isv- Taking 
advantage of this fact one can sum over the volume 5V{IVi ^ l^y in the vicinity of R in the Hamiltonian (^): 



E 

rG<5V(R) 



1 / dUr \ 1 ^ 2 



[ dt 



1 / duR\^ 1 



^^-oUr^ 



no [l-p(RPF(R), 



(5) 



E 

r'eSV{R) 



dt 



:$0Ur' 



1 fdunV 



\ dt J 



2*oUr" 



nop{Ii)SV{R), 



(6) 



1/3 

where we assume that the change in lattice constants of two end crystals is negligible, a ~ Uq , and p(R) = n(R)/no- 
Thus the Hamiltonian (m can be rewritten in the continuous medium approximation as 
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(1-p(R))Mo Uj,+(1-p(R))$ou2, 



dV. 



(7) 



The ion polarization at R is 



1 



E E '^'^r 

re<5V(R) r'eSV(R) 



SV(R) 

(1 - p(R))(7noUR + p(R)(7noVR , 



(8) 



where we assume for simplicity the same effective charges for both oscillators, and again use the fact that Isv ^ 
and consequently Ur, Vr' do not change significantly over such distances. The interaction of the polarization with the 
electromagnetic field gives rise to an additional term in the Hamiltonian: 



U,nt = P-EdV = no [(1 - p(R))uR + p(R)vr] • EdV. 



Combining Eqs. and (^) and writing out the resulting Hamilton equations for ur, vr, one obtains 

Mo ur= ^•oUr + qE, 
Ml vr= $iVr + qE. 



(9) 



(10) 



These equations of motion should be accompanied by the Maxwell equation for the electromagnetic field in the medium 

1 rf2(E + 47rP„„ + 47rPei) 



V X V X E = 



dt^ 



(11) 



where the last term describes an electron contribution to the polarization that determines the high-frequency value of 
the dielectric parameter. Another effect to be taken into account is that the local electric field on atoms that induces 
polarization is different from the macroscopic, field entering the Maxwell equations. The effective local field in a high 
symmetry crystal can be written in the form:t3L3 



E/oc — E 
Pel = noaooE/oc- 



3 ^ " 



Pel), 



(12) 
(13) 
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Substituting the effective local field into Eq. ( |lCj ) and total polarization of the volume P = Pion + Pel into Eq. (|1] 
one obtains the closed system of equations of mechanical coordinates, the polarization and the electric field. Eqs. 
(p^), (pi]), (12), and ( p^ formally resemble equations of MREI used in many papers. There are, however, important 



differences. First, the concentration parameter p(r) entering our equations is still a random function of coordinates. 
The statistics of this parameter, however, are significantly different from the original statistics of the microscopic 
distribution of impurities. The spatial inhomogeneity of the continuous function p{r) reflects the fluctuations in the 
number of impurities in the macroscopic volume SV, which are significantly reduced as compared to the fluctuations 
in the original microscopic distribution of impure atoms. More detailed consideration of the statistical properties of 
p{r) will be presented in Section VI. Second, the derivation of Eqs. (^, (pl|), (p^, and ( p^ is exphcitly based upon 
existence of a macroscopic length scale - the localization length of the local polaritons. A similar procedure cannot be 
applied to the regular phonon states, because the localization length of the local phonon states is usually microscopic. 
Finally, the requirements that the impurity-induced oscillator has its frequency within restrahlen replaces conditions 
for local phonon modes that appear in regular MREIE In this section we will assume that the fluctuation of the 
number of defects in the chosen volume is negligible and, therefore, p(R) = p. Excluding the polarization and lattice 
displacements from the system yields the following equation for the Fourier components of the electric field 

k X k X = — e(w)E„. (14) 

In this equation, e(u;) denotes the effective dielectric function of the mixed polar crystal ABi^pCp 

1 2 /, , dl dl 
e[uj)-eoo y ^, (15) 



'0 ^ ^ ^'1 ~ 

where fig ^ ~ (^"0, 1/^04)^^^ are the lattice eigenfrequencies and parameters 

£00+2 



do.i = 4^- 



SCoo Mo,! 



determine the width of the polariton gaps of the end crystals. We assume that the high frequency dielectric constant 
does not depend upon the concentration, eoo can be expressed in terms of polarizability aoo as 

= 1 



1 - ^aoo^o 



From Eq. (|l4| ) one can separate the transverse and longitudinal modes as 

fc^Ei = , (16) 

e(c.)Ejl = . (17) 

A similar expression for the dielectric function was obtained previously in Ref. |l^ within the MREI. Nobody, however, 
assumed before that the fli lies within the restrahlen of the host. It is only the concept of local polaritons that justifies 
using this dielectric function in the case when the frequency of the defect oscillator falls into the polariton band gap. 
Comparing Eq. ( [l5| ) with the respective expression of our earlier work on the 1-d model Ref. |3j, one can see that Eq. 
( p^ ) is a zero-order approximation in a series expansion in terms of the small parameter 1/Iq 1. In Section VI of 
the paper we shall demonstrate this fact implicitly. 



III. LONG-LIVING QUASISTATIONARY LOCAL POLARITON STATES IN RESTRAHLEN BAND 

The model developed in the previous section is implicitly based upon the assumption about the existence of lo- 
cal states with frequencies within the Restrahlen. In this section we shall provide additional justifications for this 
assumption. r-in 

In the original paperaH3E3 where local polaritons were introduced, it was assumed that within the restrahlen there 
is a genuine spectral gap with the phonon DOS being equal to zero over some frequency region. Such situations 
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are, indeed, possible in some crystals {GaP, ZnS, CuBr, etM). It is more common, however, that the restrahlen 
is filled with LO phonons, which are linearly coupled to local polaritons. This coupling results in the "phonon 
radiative decay" of the local polaritons (compare to radiative decay of local phonons due to coupling to light). It is 
important to emphasize, however, that the electromagnetic component of the local polaritons remains localized even 
when there is leakage through the phonon component. The phonon radiative decay broadening of local polaritons 
(and the respective life-time) is determined by the density of LO phonon states within the restrahlen. In some cases 
this DOS is large enough to suppress Jiie local polaritons. At the same time, there exists a broad range of materials 
(for instance, NaF ,NaBr, RbF, etcSIi) that have a relatively low DOS inside the restrahlen, so one could expect 
that local polaritons can survive the presence of LO phonons and affect significantly the transport properties of the 
system. In fact, the local states inside the restrahlen observed in Ref. |l| resided in the frequency region with a rather 
large DOS, and still could be observed both in reflection and Raman spectra. 

Since the fundamental assumptions incorporated into the model presented in this paper are based upon existence 
of long-living local polaritons, it is of great importance to study their life-time in realistic polar crystals. Our 
consideration of the question about the "phonon radiative decay" time of the local polaritons is based upon results 
of Ref. |2^. It was shown that the frequency of the local polariton is determined by 



6m 2 f P\\{z) + 2pj_{z) 

1 — UJ / ?; ^1 dz 



m 
6m 
m 



P||(z)-f 2px(^) ^^ 



(18) 



where m — / ™^ ; w+ and m_ are the masses of positive and negative ion sublattices respectively, 6m is the 
difference between masses of the defect and the host atom replaced by the defect, and P||(w) and p±{lli) are densities 
of states of TO and LO phonons, respectively. If the LO DOS differs from zero in the restrahlen region where we 
expect the local state to reside, the respective integral in Eq.(^8|) acquires an imaginary part equal to iTrwpn (w). 
For a light impurity {6m < 0) this equation would have a real solution cui if p\\ (lu) — O.Podolsky The presence of LO 
modes makes the solution complex valued, = lui — ij. The phonon component becomes delocalized, although the 
electromagnetic component remains localized. 

If, as we assume, j/lui 1, then lui can be considered to be the solution of Eq. ([l^ ) with the imaginary part 
dropped and 7 can be found as the first order correction to it 

where i?(wi) is a principal value of the integral entering Eq. ([l^) . The derivative of i?(wi) can be estimated using the 
explicit form of this function as uji '^"^^^'^ = 2i?(aji) + O ((wia/c)^). Finally one obtains 

7 TT / —Sm \ m_ , , , ^ 

- = T 7 — ^iPlI (^i)- (20) 

It follows from Eq. ( pO[ ) that there are three factors affecting the life-time of the polariton states. First of all, it is 
the density of LO phonon states ujip\\{u!i) in the restrahlen. As mentioned above, due to the strong dispersion of the 
LO branch, the density of states in many alkali halidesc3 between Q(^xo){0) and ^(lo){0) can be less then 10% of 
the maximum DOS at TO frequency. Unfortunately, the experimental data on the phonon DOS known to us do not 
provide its value in absolute units, and in order to obtain an estimate for this quantity in the region of interest, we 
have to rely upon some assumptions. To this end we use the Debye model for DOS of the acoustic phonons in order 
to establish a scale for the experimental results listed in Ref. This choice seems reasonable because the Debye 
model gives a fairly good description of low-temperature thermodynamic properties of crystals, and parameters of the 
model for most of the crystal of interest are established with good accuracy. Having the scale for the phonon DOS, 
one can assess the numerical value of the dimensionless quantity LUip^^lcoi). Our estimates show that this parameter 
in the spectral region of interest can take values between 0.2 and 0.6 in different systems. 

The next factor affecting the life-time of the local polaritons is the defect parameter 6m. This parameter cannot 
be assumed to be too small since we would like to have our local state farther away from the TO boundary of 
the restrahlen, and as it was shown in Refs. the frequency of the local polariton state moves deeper into the 

restrahlen with increase of the defect strength 6m. At the same time, it is clear that the factor ~6m/ {m_ + 1^+) is at 
least less then unity. The third factor m-/m+ (or m+/m— when the negative atom is replaced), can be as small as 
0.2 for RbF or even 0.06 in the case of Lil. Combining all terms we find that in a number of crystals {NaBr, NaCl, 
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RbF, etc,) the dissipation of the local polariton state is rather small 'y/uii < 0.1. In fact, it appears to be of the same 
order of magnitude as anharmonic absorption, and later in the paper we will introduce both these relaxation channels 
phenomenologically using just one parameter of relaxation. The presented estimates show that local polaritons can 
actually survive even in materials with restrahlen filled with LO phonons, and provide, therefore, a foundation for the 
model presented in the previous section. Now we can start discussing the results following from this model. 



IV. IMPURITY-INDUCED POLARITON BAND IN MIXED POLAR CRYSTALS 



A. Dispersion laws of impurity polaritons 

In this section we discuss properties of polariton excitations in our system neglecting relaxations. Effects of the 
dissipation will be incorporated in our treatment of reflection spectra in the next section. The dispersion of the 
transverse polaritons is determined by the equation 



fc = vM^-, (21) 



c 



while the longitudinal excitations obey the equation €{lo) — and are dispersionless within the present model. The 
sign of determines the structure of the spectrum. The bands of propagating electromagnetic waves coupled to 
the lattice vibration appear at frequencies where e(a;) is positive, and bandgaps arise where the dielectric function 
( |l5| ) is negative. The change of the sign of e{u}) occurs when e(w) = and when l/e(w) = 0. In the electrostatic 
approximation, the first of these conditions determines the LO frequencies lulOi while the second gives the frequencies 
of TO phonon modes, loto- With retardation taken into account for regular polaritons in pure crystals, the latter 
becomes the short wave limit of the lower transverse polariton branch, while lolo determines the k — Q frequency of 
the upper polariton transversal branch degenerate with the longitudinal phonon mode. As we shall see in this section 
for the impurity induced polaritons, the interpretation of the boundary frequencies is quite different. 
In the absence of defects, p = 0, the restrahlen polariton band stretches between transverse , 



and longitudinal. 



frequencies. Introducing defects with fl^^^^^ < Q.i < rtg^^^"^ in the system with concentrations significant enough 
to satisfy the condition I/Iq ^ 1 yet still small in the sense that p <C 1, one can rewrite Eq. ( [l5| ) in the linear in p 
approximation: 

/q(LO)2 ,,2\ ,,2 ,,2 

It can be directly seen from this expression that the impurities give rise to the band of propagating excitations inside 
the restrahl of the original crystal, with boundaries given by 

2 o2 f ldl{nl^cu^)~ldl{nl^c.^) \ 

"''^"^"H — ^¥^1 — j' 

, „ f '-^dliiu^ - nl) - ^dliuj^ ^ , ^ 

"™ ^ ' + ^ — »2 _^(^V -) ■ (23) 

The width of the band (in terms of squared frequencies) is linearly proportional to the concentration 

4eoo j2 eoc+2/Q2 _ r)2\ 
A2 - d'^iV? 3 "O 3 '^"l "O-l ^ n ^241 

A„„p-di(i2i ^^O)/ W „(TO)2X^^>"' ^^^^ 



- f^2) [^2^ _ ^(T 



provided that the impurity frequency f2i is not too close to ^^'^^ such that 
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ni"-"^^ ^-nl< -dodi^. (25) 



The lower band boundary, given by Eq. (|2^), corresponds to a dispersionless longitudinal polariton branch 

c.f°^(k)=c^,,, (26) 
while the upper one is the short wave limit of the branch of transversal excitations. An approximate dispersion law 



for these excitations can be obtained from Eq. (21) if one substitutes cu = fli everywhere except in terms containing 
impurity polariton band boundaries: 

o.(™)(k)=^,,+<5„„,^^^. (27) 
Parameter Iq in this equation is the localization length of a single local polariton with the frequency equal to lli: 

^"'^{^'"'nl-d^^i) ■ ^^^^ 

Excitations described by Eqs. (|2^) and (|27| ) demonstrate a number of peculiarities. First, one can note that 
the mutual positions of longitudinal and transverse modes are reversed compared to the regular polaritons: the 
longitudinal mode has lower frequency than the transverse one. However, if one takes the original polariton branches 
of the host crystal into consideration, the normal sequence of transverse and longitudinal modes is restored: the 
host transverse polariton branch is followed by the impurity longitudinal mode, which is followed by the impurity 
transverse mode. The last modes in the sequence are the LO mode and upper transverse polariton branch of the 
host. Second, the single transverse impurity polariton mode combines properties of lower and upper regular polariton 
modes. Indeed, at fc = this mode becomes degenerate with the impurity longitudinal mode akin to the upper branch 
of regular polaritons. At the same time the short wave limit fc — *■ of the same mode corresponds to the TO frequency 
of the electrostatic approximation, and sets the upper boundary of the propagating band similar to the regular lower 
polariton branch. 




FIG. 1. Dispersion curves of the transverse (solid) and longitudinal (dashed) optic defect modes in the restrahl for four 
different concentrations. The insert depicts the dispersion curves of the pure crystal. Frequency and wavenumber are given in 
the same units, cm~^. 
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The dispersion curves of the transverse excitations, obtained from general equation ( pT] ) for several concentrations, 
are shown in Fig. 1. The similar dispersions were observed experimentally by means of Raman spectroscopy in mixed 
crystal Ga/InP in Ref. ^ Originally, the interpretation of these observations was based upon a phenomenological 
dielectric function with multiple resonances, while it seems clear now that they provide a solid support for the concept 
of impurity polaritons. 

From the dispersion law described by Eq. ( p7|) one can obtain an expression for the group velocity of the respective 
excitations: 



which reaches its maximum value of 



c A w2A2 

^ ^ ^ imp 



c X S uj uj 

at 



(30) 



imp ' 



We introduced here a new parameter. Simp, characterizing the width of the impurity polariton band in terms of 
frequencies themselves: Simp — ^fmp/"^^- The group velocity is linear in concentration and is significantly smaller 
than the speed of light in vacuum. This is an interesting result showing that propagation of light through mixed 
crystals can be significantly (and controUably) slowed down. The smallness of the group velocity is also reflected in 
the flatness of the dispersion curves presented in Fig. 1, as well as in experimental results of Ref. 

The density of polariton states in the impurity band (IPDOS) can also be obtained from the dispersion equation, 
Eq. (H), 



1/2 

^(^) - T^Z^^^Lp—^ (31) 



27^2/3 """""P. 2 2^ 



At the low-frequency boundary of the band this DOS reproduces the singular behavior characteristics for the lower 
polariton band of pure polar crystals, while at the high frequency edge it resembles the behavior of the regular 
upper polariton band. At the center of the impurity band IPDOS can be expressed in compact form in terms of the 
localization length Iq and the linear width of the band Simp'- 

" ^^^F^- ^^^^ 

n iQUimp 

The dependence of IPDOS upon concentration is different at the edges of the band and at the center. The boundary of 
IPDOS depends linearly uponp, and tends to zero whenp decreases. At the band center IPDOS is inverse proportional 
to p, and goes to infinity when p — > 0. Such a behavior has a simple meaning - it describes the collapse of the band 
into a single local state with an infinite density. 

The linear in concentration approximation for the impurity band fails when the defect frequency Qi falls close to the 
band edge, SIq^*^'' , of the pure crystal, so that Eq. ( p5[ ) does not hold anymore. In this case, one can obtain approximate 

expressions for characteristics of the impurity band using an expansion of Eq. ( p^ ) in powers of JIq^'^'' — fli. The 
zeroth order in this parameter leads to a square root dependence of the band width upon concentration: 

AL,-^VP- (33) 

It is interesting to note that this situation was probably observed in one-two-mode mixed crystal 00/171^5.0 In that 
system at small concentration of /?i the frequencies uJn^p) and ri^^'^)(p) could not be fit with linear in p expressions, 

while uJiu{p) remained linear. This is exactly what our model predicts in the case when fli approaches Qq^'~'^: 





— SJq 


2(^0^1 r- 
3 ^' 


2 




2 £oo ,2 

32 + .^^^^' 


^{LO) 2 


= ^ - 





(34) 
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DOS in this case is given by the expression 

D{cj) 



27r2A3 ■ _ 



imp 



3/2 



/ 2 2^5/2' 



(35) 



which has the regular polariton behavior at the boundaries similar to one described by Eq. (|3^), but with a different 
pre-factor. This pre-factor significantly modifies IPDOS at the center of the band, 



31/25 



1 



(LO) 2~ 




3/2 



1/2 



^2 ^3^(LO) ^M^dl 



{LO) 



(36) 



which is now proportional to p^/^. At very small concentrations, however, this dependence will change over to the 
one described by Eq. ( |3^ ) as the condition given by Eq. (25) becomes valid again. 

The dispersion relation for the impurity TO branch becomes in this limit more complicated: 



(TO) 2., ^ 

(k) 



2Al 



3ec 



-4A4 

■imp 



(37) 



while it still has the same limits Wj-"^'^'' (k) — luu ~ fc^ for small k and ujY (k) = ujiu for large k. The crossover 
parameter, however, is no longer the localization length Iq, but rather a vacuum wavelength of light at frequency 
ri^*^: A = il^'^/c. The effective mass of the branch at the longwave boundary is much larger than in Eq. ( [27| ) 
and does not depend upon concentration. The LO frequency of the impurity polariton branch does not show any 
dispersion: 



XTO). 



(38) 



similar to the situation considered before, for k > Af^^p/ (A c^q) two branches of the impurity band run almost parallel 
showing a very small dispersion. 

The group velocity at the boundaries approaches zero at the same rate as before, but at the center of the band, the 
dependence upon the concentration is different: 




5 \ 3eoo n 



(39) 



It is much larger than the respective quantity, when the linear in p expansion is valid, and increases much faster 
~ pi/** with the concentration. 



B. Evolution of the polariton impurity band boundaries with composition parameter 

Further increase of concentration leads to more complicated dependencies of the bandwidth, DOS, etc. on concen- 
tration and other parameters. In this subsection we shall focus upon concentration dependencies of band boundaries, 
which are experimentally identified with TO and LO phonon frequencies of the electrostatic approximation. These 
dependencies were extensively studied experimentally, as we have already discussed in the Introduction, and the ob- 
jective of this subsection is to show how the concept of impurity induced polaritons provides a simple and physically 
transparent explanation for the weak mode in one-mode crystals, and for one-two-mode behavior. Analysis of poles 
and zeroes of the dielectric function ( p^ shows that the evolution of the modes with concentration is determined in 
our model by a relative position of the characteristic frequencies Jlo, ^1, ^o*^: and Of'^. The first pair 

of these frequencies are the initial phonon frequency of the host crystal, and the local polariton mode of the impurity, 
respectively. The electrostatic interaction, and the local field changes Qq to actual TO and LO host frequencies 
and Og*^. The last pair of the frequencies, fi™ and ilf*^, corresponds to TO and LO modes of the crystal made 
up of impurity atoms only. The local field, which induces the difference between initial frequencies, rio,i, and actual 
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TO, LO frequencies, is very important. In pure crystals (made of initial host atoms or initial impurities) the relation 
between f2o,i and ri™ and JIq^j 

(2 + eoo)f^L = 2f^[,™' ' + eoo^^oT (40) 

was originally derived in Ref. and its importance was stressed in Ref. The same relation exists in our model 
as well. 





FIG. 2. Four different types of tlie evolution of tlie transverse (solid lines) and longitudinal (dashed lines) optic mode 
frequencies of the mixed polar crystals with composition parameter. The leftmost and rightmost composition parameters 
correspond to pure end crystals. 



At small concentrations, p <C 1, the initial impurity related frequency fli is not renormalized by local-field correc- 
tions because this renormalization is caused by the interaction with like atoms. At small concentrations an impurity 
atom is mostly surrounded by atoms of the host crystal, and the renormalization does not occur. Similarly, for 
1 — p <^ 1, the local field turns fli into and fif while leaving Qq unchanged as a characteristic frequency of 
former host atoms. Eqs. (|2^) demonstrate that at small p the polariton impurity band inside the restrahlen of the 
host arises when J7i falls inside the host restrahlen, see Fig. 2a,c,d. For 1— p ^ 1, host and impurity atoms exchange 
their roles. In this case, should Q.q falls in between Vl'"^^^ and it gives rise to the impurity polariton band 

induced now by the "host" atoms. Hence, if both f^o"^*^'' < ili < VL^^'^^ and ri^^*^-* < i7o < ri^^*^^ conditions are 
satisfied at the same time, then the polariton impurity bands exist in both p — > and p ^ \ limits. This situation is 
shown in Fig. 2a. This is one-mode behavior with a "weak" mode. One observes strong TO and LO modes of the 
original crystal smoothly evolving into those of the end crystal, while inside the restrahlen a weak additional TO and 
LO impurity polariton modes arise with vanishing TO-LO splitting at both ends of the impurity range. This peculiar 
behavior occurs due to the interplay between the external electric field and the polarization, affecting the local field 
(|l^). Experimentally, this type might be realized in the materials where the polariton gap is sufficiently wide, which 
may be the case for many alkali halides. However, because of the weakness of these modes they are vulnerable to 
any kind of dissipation, as we shall see in the next section. This fact can explain the absence of the defect mode 
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in "classical" (no weak mode inside the restrahlen) one-mode mixed crystals. At the same time, we can relate the 
features presented in Fig. 2a to the weak mode observed in the Ba/SrF2, Ca/SrF2 and some other one-mode crystals 
(see Discussion). Previously, in order to reproduce the weak feature observed in spectra of these crystals, one had 
to use models with tens of fitting parameters (see Refs. |l|,^ and references therein). The concept of the polariton 
impurity band, presented here, gives a transparent and quite general explanation of this type of spectra. 

If the polariton impurity band exists only in one limit, for example, f^o^'^'' < ili < r^o^*^' for small p, but flo lies 
outside the interval (fi^^'^-', fi^^*^-*) whenp ~ 1, our model predicts the onc-two-mode behavior (Fig. 2c, d). Two types 

of one-two-mode mixed crystals appear depending upon whether Qq falls below the interval (il^^'^\ (Fig. 2c), 

or above it (Fig. 2d). In the first case, the lower mode weakens with the concentration, while in the second case, 
the upper one does. At p ~ the situation remains qualitatively the same as in the one-mode case. At p — > 1, the 
splitting between new modes does not vanish, and they form the new restrahlen of the pure crystal ai p — 1. Strictly 
speaking, in this limit we cannot justify our model, since there are no local modes at this end of the concentration 
range. Quantitatively however, we can still explain the one-two-mode behavior as a result of the presence of the local 
polaritons at small p, and the absence of them at p close to one. If in neither limit the impurity polariton band appears 
inside the gap {fli <f. [^^'^^ and fio ^ (fij'^'^^ fij^*^^)), our approach is not applicable at any concentration. 

Although the absence of the impurity bands at end concentrations, in our opinion, suggests two-mode behavior (Fig. 
2d). 



V. TRANSMITTANCE AND REFLECTANCE OF EM WAVES FROM MIXED POLAR CRYSTALS 

Having obtained the effective dielectric function of the mixed polar crystal, one can study their reflection and 
transmission spectra for a number of different geometries. Among them, the normal reflection spectrum from a semi- 
infinite sample and normal transmission and reflection from a slab of finite dimensions are of particular interestpsince 
they are a primary source of information about optical properties of crystals studied in numerous experiments.EJa 



A. Reflection spectrum of semi-inflnite mixed polar crystals 



In this section we phenomenologically incorporate damping in the dielectric function (|15|) substituting ii^^i'^ for lo in 
the resonance terms of Eq. (|l5|) . As we already discussed it is sufficient for our qualitative purposes to assume that the 
damping parameter is the same for both host and impurity related modes. We can note, however, that as numerical 
calculations demonstrated, properties of each mode are determined primarily by its own relaxation parameters, and 
effects due to damping in the other subsystem are usually negligible. This means that, in principle, it is possible 
to experimentally determine relaxation parameters for each of the participating oscillators independently. For our 
calculations we chose the value of 7 such that 7/rio ~ 0.1, which is a rather conservative estimate. For typical polar 
crystals the relaxation parameter ranges from less than 0.01 to 0.1, and since the relaxation in the impurity 

subsystem due to coupling to LO phonons inside the restrahlen was estimated in Section III as less than 0.1, our 
choice for this parameter seems quite reasonable. 

After accounting for the dumping, the effective dielectric function (p^) becomes 



e(^) 



2 

1+3 



■ up- -I- 2i7cj 



d\ 



d^ 
"0 



di 



(41) 



For normal incidence, the reflection coefficient from the semi infinite crystal is equal to 



.1/2 



(a;) 



l + ei/2(a,) 



(42) 



In the absence of damping, the original dielectric function ( |l5| ) has a peculiar property, which is specific only for the 
model with the impurity polariton band within the restrahlen of the host. Inside this band, the dielectric function goes 
from zero to infinity, and hence, necessarily passes through unity. At a frequency where this happens, the reflection 
coefficient must become zero, since for this frequency the medium becomes transparent. At small concentrations, this 
frequency, lot=i-, is determined by the equation: 
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and since ^ Ai = fii/c it lies slightly off the center of the band, closer to uju. When relaxation is accounted for, the 
zero of reflection is not reached, but the reflection still can have a minimum at a certain frequency. The magnitude of 
reflection at this frequency is determined by relaxation, and can be used for independent measurements of the latter. 

We used typical values of the parameters in Eq. (|4^) to plot normal reflection spectra for semi-infinite mixed 
crystals. The optic frequencies J^q™'^*^^ of such systems are about a few hundred cm^^, the widths of their polariton 
gaps range from ~ 10% for III/V group polar crystals, up to ~ 30% for alkali halides. The high limit dielectric 
constant lies within the range of ~ 3 — 5. Fig. 3 presents three graphs corresponding to three different types of 
spectra, which can be described within the model of impurity polaritons. Each graph shows curves obtained for 
p = 0%, 25%, 50%, 75%, 100%. 




FIG. 3. Normal reflection spectra from semi infinite mixed crystals. Three graphs correspond to one-mode (graph (a) in Fig. 
2) and two types of one-two-mode (graphs (c) and (d) in Fig. 2) behaviors. The parameters used to generate these graphs 
(including absorption) are typical parameters for polar crystals. 

Fig. 3a depicts the reflection spectrum of the one-mode crystal. As it can be seen from the plot, for all concentrations 
there is one dominating absorption band. Nevertheless, for intermediate concentrations one can notice a weak mode 
inside this band. This mode is quite weak in accord with the discussion of the previous section, and can be smoothed 
away by the absorption. Whether this mode will be observed in a concrete material depends upon the interplay of 
several parameters. At the same time, since we have used realistic values of parameters characteristic for one-mode 
group of mixed crystals, our calculations show that the impurity polaritons can indeed be used to explain this feature 
of one-mode type spectra. Reflection spectra corresponding to one-two-mode behavior are shown in Figs. 3b and 3c. 
There are two types of such spectra, which are very much alike. At small p, the spectra look akin to the one-mode 



14 



type, but with an increase of the concentration two modes appear, with one-mode growing stronger and the second 
one diminishing. Spectra 3b and 3c correspond to phase diagrams shown in Figs. 2c and 2d, respectively. It is seen 
that the presence of damping does not prevent one-two-mode behavior predicted by our model to be observed in 
reflectance experiments. 



B. Slab of finite dimensions 



In this subsection we consider normal reflection and transmission spectra of a mixed crystal slab. The width of the 
slab is assumed to be much greater than the average distance between the impurities, so that our averaging procedure 
can be applied. At the same time we do not consider samples thicker than several wavelength, so that the damping 
does not suppress transmission completely. 

The transmission coefficient through a slab of width L for normal incidence is given by 



T(c.) 



1/2 



l)2cos(^ei/2L) 



2ei/2 cos (^e^/^i) _ ^(g + i) sin (^e^/^L) 



(44) 



In the absence of damping, the transmission coefficient turns to unity when either e = 1 or we^/^L/c = 7r(m + 1/2), 
where m is an integer. The flrst case corresponds to the frequency determined by Eq. (^), when the medium becomes 
optically transparent, while the second condition corresponds to the usual geometrical resonances. 
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FIG. 4. Transmittance spectrum through a slab of one wavelength width with typical parameters, plotted as a function of 
frequency. The lower graph shows the full frequency range, while the upper one is restricted to the polariton defect band at 
p = 10%. 
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The width of geometrical resonances decreases as l/e(w) with lo Wi„ because of the divergence of the dielectric 
function at this frequency. Fig. 4 depicts the transmission spectrum of the 10% impure slab of the width L = A in the 
absence of absorption, where the lower plot shows the broader frequency interval covering the entire restrahlen of the 
host crystal. One can see the impurity induced polariton band inside the forbidden gap at the frequency ^ 180 cm~^. 
The wide flat-top resonance at 181.7 cm~^ corresponds to the frequency (^3|), where the dielectric function becomes 
one, and all the other peaks represent geometrical resonances. 




100 200 300 400 

a 

FIG. 5. Transmission (solid) and reflection (dashed) coefficients of a thin mixed crystal at p = 20% in the presence of 
absorption. 



Different natures of the resonances affect their response to damping. The narrow comb of geometrical resonances 
washes out first even when damping is relatively small. This occurs due to the fact that the geometrical resonances 
are in essence standing waves that experience multiple reflections and therefore are strongly damped. The e = 1 
resonance, to the contrary, is much less affected by the relaxation because polaritons pass through the sample only 
once. Fig. 5 depicts the transmission (solid line) and reflection (dashed line) coefficients for a narrow slab with 20% 
impurities and the same set of parameters, which was used to generate Fig. 3c. It is clearly seen, that even despite 
of high absorption rate j/^o — 0.1, the peak in the transmission coefficient survives, and can be associated with 
the minimum in reffection. It allows us to suggest that reflection anti-resonance observed in many systems can be 
associated not only with the impurity induced absorption, but also with the impurity induced transmission. A similar 
transmission maximum was observed experimentally in CuClS^ where the role of defects was played by some Cu 
atoms occupying off-center positions. The maximum in the transmission is direct evidence of the impurity-induced 
polariton band. It would be of great interest to carry out transmission measurements for different groups of mixed 
crystals in order to verify predictions of the current paper. We believe that time-resolved measurements could also 
provide important information, particularly regarding group velocities of these excitations. 



VI. SCATTERING IN THE IMPURITY-INDUCED POLARITON BAND 



In previous sections we considered optical properties of mixed polar crystals neglecting ffuctuations of concentration 
function p{r). Results obtained for the one-dimensional model considered in our previous paper, Ref. |4[ suggests 
that this approximation is a zero-order term in the series expansion in terms of the parameter 1/Iq <C 1, where I is an 
average distance between defects and Iq is the localization length of the single local polariton. In this section we take 
into account ffuctuations of concentration in the first non-vanishing order in terms of the parameter Sp — p(r) — p, 
and show that the actual small parameter of this approximation is indeed I/Iq <^ 1. The main result of this section is 
to demonstrate that the scattering of impurity polaritons induced by concentration fiuctuations is decreasing with an 
increase of the average concentration, and is actually negligible on the background of absorption. This fact provides 
a justification for the results obtained in the previous sections of the paper, where these fluctuations are neglected. 

We shall employ the Green's function formalism in order to calculate the scattering length of the impurity-induced 
polaritons due to concentration fluctuations. The Green's function of the inhomogeneous dielectric medium can be 
presented as 
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SafsSir - r') 



(45) 



where Gaj is the Green's function of the system with homogeneous dielectric function e{uj,p) given by Eq. (|l5|). This 
zeroth order Green's function can be written down in fc-space in terms of the projection operators e^^ = — kakp/k^ 

and e||,g ~ kakp/k^ as 



1,2 _ jU2 



(46) 



where fco(w) = w^e(w)/c^. Since we search for the leading corrections to the Greens' function of the system, we expand 
e(a;,}3(r)) in terms of 5p{y) and keep only the linear term 



£('^)-P(r)) - e(^,p) = 



uj"^ 5e{uj,p + 5p{v)) 



5p{v) 



Spiv) 



5p(r)=0 



= i^^i^.p) ■ Sp{r). 

Statistical properties of the random function Sp(r) can be described in this approximation by its momenta up to the 
second order: 



{Sp{r)) = , 
{Sp{r)Sp{r')) =K{\r-r'\). 



(47) 



The second order correlator depends only upon the distance between two points in space, since the system is assumed 
to be homogeneous and isotropic on average. After standard transformations to fc-representation one obtains 



G;^(k) 



d^k' ( 



(48) 
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1^12 _ 



S{k-k') 



where S'(k) is the Fourier transform of the correlator 
and longitudinal mass operators 'E'^{k), I]ll(fc): 



- fc2 - J:^{k))eip - [kl + Y.\\{k))eip, 

\. The new Green's function can be expressed using transverse 



Ga^(k) 



k"^ -kl^Y.^[k) fc2 + sll(fc)' 



(49) 



Real parts of the mass operators determine the renormalization of the spectrum, and are neglected below. We shall 
only consider imaginary parts of the transverse mass operators, which determine the scattering length (mean-free-path) 
of the transverse modes, and in the lowest approximation read as 



Im E^(fco(a)),cj) i — ko 



S ( 2fco sin - ) c?cos 9, 



(50) 



where 9 is scattering angle between k and k'. The explicit form of the integral scattering cross section depends on 
the particular choice of the correlator (^). In our case however, when the wavelength of the considered excitations 
is assumed to be much greater than the characteristic size of the inhomogeneities, the difference between different 
choices of the correlation function is mostly reduced tOj-a. numerical factor of the order of unity. One can choose the 
correlator, for example, in the standard Gaussian formO 



K{\v-v'\)^{6p') e-|--^T/'' 



(51) 



with its respective Fourier transform 
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S{k) 



(52) 



In the spirit of our general approach, the concentration fluctuations must be considered with regard to the smoothing 
volume SV, so that the respective variance (Sp^) is calculated using the Poisson distribution of independent impurities 
inside the volume SV: 



N{SV) 



(53) 



The correlation length Ic, then should be identified with the Z^y-smoothing length we employed in Eqs. (||) and 
(^. However, as will be seen below, the same results can be obtained if one considers initial distribution of discrete 
impurities and chooses the interatomic distance a as the correlation length. In any case, fegZc ^ 1, and the expression 
for the complex wave number fc(w) determined from the pole of the respective Green's function takes the form: 



Using Eq. (m) and the identification of 1^ one can obtain for the scattering length: 



i: 
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uj^ Je(w,p(r)) 



n 2 



(5p(r) 



(54) 



(55) 



In the limit of small concentrations, p <^ \, the derivative of e{uj,p) with respect to p can be evaluated in different 
regions of the spectrum. 

To assess the value of the scattering length inside the impurity band we pick the center of the defect band oj^ — 
{ufi + Lo1^)/2. For small concentration one obtains 



K' i^c) 



where 



ai 



ttV^ 2 

4 i^,{n,f' 
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P 



(56) 



4e, 



3(2 + eoo) " 2 + e 



does not depend upon the concentration. In the last expression the anticipated small parameter l/lo has appeared 
raised to the third power. Taking into account typical values of this parameter for realistic mixed crystals, one 
can see that this scattering is completely negligible. One may also note that this scattering length increases with 
concentration, which has a simple physical explanation - the greater the concentration the greater the overlap of 
individual local states, and the closer the system is to a uniform continuous medium. 
For frequencies from host bands, the situation is qualitatively different: 



where 




(57) 
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and 
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14 {n^^^'^'-n^o^''^') 

(niLO) 2 _ fl{TO) 2) (fliLO) 2 _ ^2 j 

are constants. This scattering length is much smaller then for impurity polaritons, and is decreasing with the increase 
of concentration, which is quite a regular behavior. Indeed, the perfectly ordered at p = system experiences 
increasing scattering due to the impurities. This difference between scattering properties of regular and impurity 
polaritons explains why our approach can be justified to describe properties of the former, while the latter requires 
more elaborate treatment of disorder. 



VII. DISCUSSION 



In this paper we suggested that certain features of optical spectra of mixed crystals in the restrahlen region can be 
explained if, along with standard optical phonon vibrations, one introduces additional impurity induced oscillators 
with the frequency falling into the restrahlen of the host crystal. Justification for this assumption -cpmes from the 
concept of local polaritons introduced earlier in Refs. ^J-^ and studies of one-dimensional modelsJ^Jo which showed 
how local polaritons develop into an impurity induced band. Local polaritons are significantly different from local 
phonons, which were extensively studied in connection with optical properties of mixed crystals. First, they are 
assumed to exist in the restrahlen region, which in many cases is filled with phonons. Therefore, it is not exactly a 
spectral gap, which is necessary for local states to exist. However, the dielectric function within the restrahlen remains 
negative, which means that electromagnetic excitations cannot exist in this frequency region, which is therefore a 
gap for electromagnetic excitations. This situation is reversed compared to the case of pure phonon gaps, which 
are deprived of phonon states, but contain electromagnetic ones. Because of non-zero electromagnetic DOS, local 
phonons acquire their electromagnetic radiation width, and because of non-zero phonon DOS, the local polaritons 
acquire their phonon "radiation" width. A significant difference between these two situations is that the density of 
electromagnetic states in the phonon gap is so small that the radiative broadening of local phonons is negligible, while 
the density of phonon states in the restrahlen is large, and local polaritons may or may not survive it. Analyzing 
knowp_|phonon DOS, we found that for some crystals like GaP, ZnS, CuBr, ZnTe, Cul, SrF2, BaF2, -P6F2, UF2, 
CaF2L]i there are actually frequency regions with zero phonon DOS, though often quite narrow ones. Therefore, the 
right impurity could in principle give rise to a local state considered in Refs. M^7. However, restrahlen region of a 
much broader class of crystals like NaCl, NaBr, KCl, RbCl and many otherscll are filled with LO phonons whose 
DOS within certain regions is relatively small. In the present paper we studied the life-time of local polaritons due to 
interaction with these LO phonons, and found that under regular circumstances this life-time is no shorter than the 
one due to anharmonicity. We argued, therefore, that local polaritons can actually survive interaction with restrahlen 
phonon states, and contribute to the optical properties of the crystals. Moreover, the presence of some local states 
within the restrahlen was confirmed experimentally in Ref. where neutral dopants {S, Sn, Te) give rise to the local 
electron-phonon state at the frequency slightly below LO of GaP. What is interesting is that the region where these 
states reside, has a relatively high density of phonon states, which obviously did not preclude them from existence. 
All these arguments justify the use of local polaritons to describe properties of the restrahlen of mixed crystals. 

The second important property of local polaritons is that their spatial extent is of the order of optic wavelengths 
in the restrahl (~ 10^'^ cm), which is much larger than the size of local phonons (several interatomic distances). 
This fact means that even at residual concentrations of impurities, local polaritons overlap, forming a well developed 
band. We showed in this paper that this band can be described using the continuous medium approximation for the 
impurity subsystem. Neglecting fluctuations of impurity concentration, we derived an effective dielectric function for 
our model, and used it to analyze the structure of optical spectra of mixed crystals. This dielectric function describes 
new impurity induced polariton bands, which arise inside the restrahlen of the host crystal. 

The first problem we set out to consider using tbapConcept of impurity polaritons was weak features in the spectra 
of so called one-mode crystals BaJ£rF2, Ga/SxF2,iM3 ZnCdS, Mn/ZnTep and the one-two-mode beliavior of a dit. 
ferent group of crystals In/GaP^ PbSe/Te}d KmbljBB RbBr/Gl, GaAs/Sb, InAs/Sb, AgBr/Gl^ In/GaAsH 
The existing descriptions of these types of spectruma require a great number of fitting parameters. Our model allowed 
us to explain these types of spectra naturally as manifestations of the impurity polariton mode, which reveals itself 
differently depending upon the relation between fundamental frequencies of crystals at both ends of the concentration 
range, and the frequency of the local polariton. Introducing relaxation in a phenomenological way, we considered 
reflection and transmission spectra, and demonstrated that our model survives rather strong damping, and reproduces 
spectra closed to experimental observations. For rather thin samples our model predicts a resonant enhancement in 
the transmission at the frequency of the impurity polaritons, which accompanies an anti-resonance in reflection. The 
latter was observed in many papers,™3 but was mostly attributed to the impurity induced absorption. The only 
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transmission measurements known to us were performed on pure CmC/,e3 where transmission was found to exhibit 
the maximum similar to one predicted in our paper. CuCl is a peculiar materiaL-since Cu atoms at low temperature 
can occupy several non-equivalent positions, thereby creating internal defects.E^S These off-center Cu atoms can be 
responsible for the impurity polariton band, and therefore, transmission spectra of CuCl can be considered as the 
first evidence of this band. 

Since we took the retardation into account, we were able consider not only boundaries of the spectra but also 
dispersion laws, DOS, and group velocities of the impurity induced polaritons. One of the most remarkable properties 
of these excitations is that their group velocity is proportional to the concentration, and can be thousands of times 
smaller than the speed of light in vacuum. The smallness of the group velocity makes dispersion curves of the 
excitation look almost flat. Rather similar dispersion curves were measured experimentally in Ref. 28 in one-two- 
mode mixed crystal Gao ^oIriQ ^oP with the use of Raman spectroscopy. It would be interesting to carry out additional 
steady state and time-resolved experiments in this material, which could verify predictions of our theory and provide 
more solid support for our concept of impurity-induced polaritons. Our approach also allowed us to study scattering 
of impurity polaritons due to fluctuations of concentrations. We found that their scattering length is significantly 
different from the similar characteristics of regular polaritons of the host material. The scattering length of impurity 
polaritons is proportional to the concentration, making the scattering less efficient with an increase of concentration. 
Quantitatively, the scattering length is very large, much larger than the attenuation length due to inelastic damping, 
making scattering due to concentration fluctuations negligible for impurity polaritons. This finding is in agreement 
with the results obtained for the one-dimensional model,0 and provides a firm foundation for our approach. The 
scattering length for regular host polaritons, at the same time, is inversely proportional to the concentration, and is 
rather short. Therefore, studying optical properties related to these excitations requires more elaborate theoretical 
approaches used in many papers on the subject (see, for instance, Ref. pO| ). 

We hope that the present results will revive an interest of experimentalists in the properties of restrahlen of mixed 
crystals. In our opinion, it would be interesting to study transmission spectra through thin slabs of one-two mode 
crystals in both sw and time resolved experiments. Such experiments could provide additional insight into properties 
of impurity polaritons, and elucidate their dynamic properties, such as group velocities. Comparing experiment and 
present theoretical results, one can obtain additional information about the material parameters of these systems. 

Note added in proof A.J. Sievers drew our attention to an alternative explanation of the IR reflection spectra of 
CuCl. In Refs. [39], an additional absorption line inside the gap was attributed to the interaction of TO phonons 
with acoustic phonons arise from by strong anharmonic terms in the potential of Cu atoms, rather than due to an 
additional polariton band associated with the off-center ions. Although there exists strong evidence in favor of the 



model discussed in Refs. 40, it does not explain, however, the enhanced transmission observed in Ref. pC 
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